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; Abstract 

An averaged generating function for coloured hard-dimers is being inves- 
tigated by proving estimates for the latter. Furthermore, two different enu- 
merating problems and their distributions are studied numerically. 

in 

m ■ 1 Introduction 

cn 

' Coloured hard-dimers have recently appeared in the context of causally trian- 

gulated (2+l)-dimensional gravity. Using special triangulations of spacetime, 
it was shown in [T] that the generating function of the one step propagator can 
as well be expressed through the generating functions of coloured hard-dimer 
configurations. 

In this article we go in a slightly different direction and study the behaviour 
of an averaged generating function for coloured hard-dimers. In this manner 
we touch some combinatorial problems, which so far have not been tackled. 

The paper is organized as follows. In section 2, we define the notion of 
coloured hard-dimers, give the statement of the problem and prove estimates 
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for the averaged generating function. In section 3 additional material for 
numerics around coloured hard-dimers is presented. 

2 Description of the model and results 

Let £/v be a sequence of length iV of blue and red sites on the one-dimensional 
lattice Z. Let us define the notion of a coloured hard-dimer on £w. a coloured 
hard-dimer is a sequence of red or blue dimers, which cannot intersect each 
other ("hardness" property). A dimer in turn is an edge linking two nearest 
sites of the same colour. In FigQ] an example of a coloured hard-dimer is 
given. 



Figure 1: A hard-dimer, A = 13 

For any configuration the generating function associated with it is of 
the form 

Z/. N (u,v,w) =J2u nb ^v nr ^w nbAD \ u,v,w£ (0,oo), (2.1) 

D 

where D is an hard-dimer on n b {D) and n r (D) denote the number of 
blue and red dimers respectively on D, and n br (D) the number of crossings 
between dimers and sites of different colour. Moreover, let us denote by %(D) 
and 7r(-D) the number of blue and red sites of D respectively, which do not 
belong to any dimer. We call such kind of sites "single points". Naturally the 
following constraint 

2n b (D) + 2n r (D) + n br (L>) + lb {D) + lr {D) = N (2.2) 

holds. It is useful to define further variables t and s, where t corresponds to 
the number of sites occupied by dimers, and s to the number of dimers: 

t = N - lb {D)- lr {D), s = n b (D)+n r (D) (2.3) 

In the above example n b (D) = 1, n r (D) = 2, n br {D) = 3, 7&(-D) = 3 and 
lr(D) = 1. 
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For any N, we want to prove estimates from above and from below, for 
the mean of the generating functions Z^ N defined in (|2.ip . In achieving this 
we calculate explicitly, for a given N, the sum of all the generating functions 
Z^ N over £jvj rescaled by the factor 2^, i.e. the number of all the sequences 
£tv. We get so an "averaged" generating function 

2JV zZ Z Zn( u > v > w )- 

The main result of this section is the following: 
Theorem 2.1 For any N, the following estimates 



Ci C /A + 1\ N 1 ™ 1 1 xiV 

CAT 

(2-4) 

hold, where A = \/u + v + and C\,C are constants not depending on N, 
given as 

Ci = e 6*/—, (7 = — 1 A 



2 N VN 



vr' 2 V Vu~+^+f + l 

Remark 2.1 iVofe that A > 0, so that the upper bound in \2.1$ is always 
exponentially divergent, as N — > oo, whereas the lower bound increases expo- 
nentially for A > 1 and large N, i.e. for the class of parameters u,v,w, such 
that \Ju + v + j > 1, otherwise it tends to 1. 

In order to prove the main result, we first need a Lemma, where we get 
an explicit formula for the average jjq- Y1^ N , for any N. We are able to 
get an exact formula by summing over the variables t and s defined above in 
L31). 



Lemma 2.1 The average jw^2^ N Zg N has the following explicit expression, 
for any N : 

where [■] denotes the integer part. 
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Proof of Lemma 2.1: Consider a hard-dimer D of length N, and set 
n b (D) = k, n r (D) = h, n br (D) = m, 7b(D) = 7& and 7 r (D) = 7 r . By 
(|2.2p and (|2.3p the following equalities 

2A; + 2ft + m + j b + j r = N, t = 2k + 2h + m (2.6) 

hold. 

First we fix the number of blue and red dimers and that one of single points 
7b and j r , without assigning any length to dimers, even though, by (|2.3p . the 
variable t on D is done. In this manner the dimers of the same colour are 
considered indistinguishable. Then we calculate all possible permutations of 
k blue dimers, ft red dimers, 7& blue and 7 r red single points, i.e. 

(fc + /i + 7fe + 7 r )! 
k\ h\ 7ft!7 r ! 

Now, one has to distinguish two cases: k = h = and k + h ^ 0. If k = 
h = 0, then m = 0, because this means that there are no dimers. Therefore 
lb + 7r = N and the combinatorics is 



N 'AT 



V 



E = 2 

In the other case (k + h ^ 0) we assign, for any given permutation of 
coloured dimers and single points, all the admissible lengths to each dimer, 
taking into account the fact that t is given. The number of all these combi- 
nations is then 

f m + k + h — 1 N 
k + h-l 



Therefore we have 



N 



E E 



^2z^ N (u,v,w) 

Civ 



(fc + ft + 76 + 7 r )! fm + k + h - 1 



fc! hi 7 b ! 7r ! V k + h-l 

k,h,m: lb + lr =N-2k-2h-m 10 ' v 

2k + 2h + m=l 



u k v h w m 



(2-7) 

From (|2,6p we have that A; + + 7^ + 7 r = N — k — h — m and = N — 2k — 
2ft — m — 7b and hence 

^Z Cjv (u,v,u;) = 
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N 2k 2h m (iV-fc-fe-m)! fm + k + h-l\ k h m 

^ ^ k\ h\ lh \(N-2k-2h-m- lb )\ I k + h-1 ) 

k,h,m: 
2fc + 2/i+m = l 

^ k\ hUN - 2k - 2h - m)!V + / 

2fc+2h + m = l 

N-2k-2h-m , AT , , v 

^ f * - 2k - 2h - m ) (2.8) 

In (|2.8p we have multiplied and divided the generic term of the sum by (N — 
2k — 2h — m)\, so that the binomial formula appears 

N -2k-2h-m /JV — 2k — 2h — m\ = ^2k-7hr-m 
and ^2^ N Z^ N becomes 



2 N 



E X (N — k — h — m)\ (m + k + /i — 1\ /u\ fc /f \ h fw 

k\h\(N - 2k - 2h - m)\ \ k + h-1 ) \l) VlJ V"2 



fe ,h,m 
2k + 2h + m=l 



(2-9) 

Note that the scaling factor 2^ naturally appears in (|2.9p . Now, by mul- 
tiplying and dividing the generic term of the sum in (|2.9|) by (A; + h)\ we 
get ^ 

^ ^ ( 'N — k — h — m\ (m + k + h — 1\ (k + h\ fu\ k {v\ h {w\ m 



k,h,rn: 
2k + 2h+m = l 



k + h \ k + h-1 J V k 



We perform the variable changements s = k + h and t = 2k + 2h + m, where 
s and t were already defined. We get thus 



N l 2 J 



t=l s =l v 7 v 7 fc=0 V 7 



1 +£|rriC;:;0(T)*(!r <-> 

In the last sum in (|2.10p we used again the binomial formula 



y% i s \ ( ii \ >■ ( r \ f a + r 



k=0 



k U 



Therefore we get formula (|2,5p where only the indices s and t appear. The 
Lemma is so proved. □ 

Proof of Theorem 2.1: First we write the sum on the right hand side of 
formula (|2,10p in a more suitable form: 

N-t + s\ ft-s- l\ s( 2s ) fN\fN-2s\ 



s 



s-1 J N(, N -\) \2sJ V t-2a 



(2.11) 



In fact, it holds 

'N - t + s\ ft - s - 1 



s J \ s — 1 

{N-t + s)\ (t-s-l)\ N\ (AT -2s)! (2s)! 
s\(N -t)l ' (s - l)\(t - 2s)\ ' W. ' (N - 2s)! ' (2s)! 



s!(s-l)! 



iV\ fN - 2s 



A[! V2s/V t — 2s 

*( 2 ;) w ^-2« n 

Now we estimate the binomial coefficients ( 2 S S ) and (^^-J from above and 
from below, more precisely, first we prove that 

^ < iffl < 2 ' 2 ' ( 212) 

where C\ = e~e y^- In order to estimate the binomial coefficient ( 2 g s ) we use 
Stirling's formula 



s! = s s e" s \/27rs e As 



where 
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We may write 

2 ' S ^ 2 e A 2s~ 2A s (2.13) 



s J ^/^^s 

Taking into account the following inequalities for the errors 

-7T < ~ 2X S < 

6s 

we get estimates from above and from below for s( 2 s s ) 

' " 2 2s <s( 2S ) <2 2s v^ (2.14) 



7T 



For the binomial coefficient (+^._\) we use the fact that its maximum is 
assumed at t — s — 1 = [^f^] ■ Hence 

/ N-l \ fN-1 

We can assume for simplicity that iV is even. The case of iV odd can be 
treated with little more effort. Then, as before, we use Stirling's formula 

- 1\ N\ 

N 



f-V (f-l)!(f)! 
N N\ < N2^_ = 2N IN 



2 (f )|2 - 2 Jn V 2 



2 

Thus we obtain 

that yields with (|2.14p the formula (|2.12p . Therefore we are able to estimate 
the mean J2^ N from below and above using the same series 

2 N y/N ~i \2s J \ t — 2s J ' \2 



t=l s=l 
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1 



N I5J 

1 + — >:>: 



N\ ( N - 2s 



2s 



t=l s=l 

It remains to evaluate the series 



t - 2s 



(u + vf 



t-2s 



N I2J 

EE 



iV\ / AT - 2s 



2s 



t=l s=l 

We then state the next Lemma 



t- 2s 



(u + v) 1 



Lemma 2.2 T/ie following equality 



w 



t-2s 



(2.15) 



(2.16) 



A? 



EE 

t=l s=l 



N\ fN - 2s 



2s 



t - 2s 



{u + v) 



w \ t—2s 

~2 



( 1 w \ n ( , w 

\^fu + v+- + lj + (-yju + v + - + l 

holds, for any N . 



N 



(2.17) 



Proof of Lemma 2.2: Consider the sum (|2.16[) for s and t starting from 
and denote it by *. Note that in * there is an extra term at t = s = which 
equals 1. In this case we can apply the following combinatorial equality (see 

0) 



E 



n \ l n — v 



v ) \k — v 



(t + 1) 



(2.18) 



where the sum is over any u such that < v < k < n. The summation with 
respect to s in * coincides with the sum (12. 181) restricted to the even indices 
v. Owing to this fact we rewrite the sum over s in a more suitable form: 



E 

s=0 



N\ fN - 2s 



2s 



t - 2s 



E 



1 + {-iy (n 



1 * 
-Y 

2 ^ \ s 

s=0 



N -s 
t - s 



(\/u + vY 



N\ f N — s 
t - s 



2^/uTv 



w 



+ 



1 * 

-Y 



N\fN-s\/ 2^/U^\ /w\t 
t-s { w ) V2. 



We apply now formula ()2. 18j) 



=0 



2s fW\ t - 2s 

2. 



1 (N 

2 V t 



+ 1 + 



2y/u+~ 



+ 1 



1 (N 

2 V t 



(2.19) 



Summing the last term in f|2. 19j> over t we get again a binomial formula 



1 N 



t=Q 



w\t ( , u>\* 

yju + v + —j +y-\/u + v + — 



( i w A N ( i w -, xAr 

(y/u + v + j + lj + [~yju + v + - + 1 

The Lemma 2.2. is so proved. □ 

Now we are able to finish the proof of the Theorem 
Upper bound: Taking into account (12.150 and (12.170 we obtain an upper 



bound for ^ % 



1 1 / . W \N 

1 - -7= + —F= [\ / U + V + - + 1 

y/N 2^/N V 2 



1 + 



AT 



Since 



< 1, 



v 7 ^ + v + f + 1 
we have proved the upper estimate in (|2.4p . 

Lower bound: Analogously, we obtain as lower bound for Egjv 



Ci Ci 



, U> \ N f i w \ N ' 

+ u + - + lj + (-Vu + v + - + l 
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Ci + Ci 



2^x7^ 2VN 



y/u + V W 1 

^ + 4 + 2 



N 



+ 



Ci Ci / V^+w to 1 
+ 2~7^ V 2 + I + 2 



1 + 



+ f to 1 
^ + 4 + 2 



-V^+u + t + 1 



A'l 



V^+f + 1 



N~ 



Therefore we get the lower bound (|2.4p . with the constant 



(2.20) 



C = y for -^IT^ + ^ + 1 > 



C=^|l + 



-VnT^+ f + 1 



Ci 



w + 2 



2 VV^+^+f + 1/ ' 



for — + ^ + 1 < 0. The Theorem is so proved. 
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3 Numerics for coloured hard-dimers 

In this section we provide some supplementary material regarding hard-dimer 
combinatorics. We consider two enumeration problems which are different 
from that in section 2. Numerical techniques are being employed and the 
results may indicate interesting directions for further study. In the first part 
we count the number of hard-dimers for each configuration £jv> N fixed. The 
so obtained statistics is then shown in histograms. 

In order to explain the algorithm we consider a particular configuration as 
shown in Figj2j 



Dimers for a partic- 
ular configuration of 
sites are indicated as 
closed intervals. 



Figure 2: Dimers 

Recall that a hard-dimer on a sequence of red and blue sites is a sequence 
of dimers which are non-overlapping. In order to count the number of hard- 
dimers on the fixed sequence £at of red and blue sites, one can proceed as 
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follows: 



Step 1. Generate a list of all dimers, i.e. codify their colour as well as their 
starting and end points. E.g. for the dimers in Figj2] this is achieved by 
defining the set 



{{1, 4, r}, {2, 3, b}, {3, 5, &}, {4, 6, r}, {5, 9, b}, {6, 7, r}, {7, 8, r}, {8, 10, r}} 



Step 2. For each I £ {1, . . . , |-E-D|} let K(l) be the set of subsets of ED con- 
taining I elements. Do the following for I = 1 to I = \ED\. Scan the elements 
of K(l). If for an element in K(l) the dimers are non-overlapping, increase 
the number of hard-dimers by one, otherwise keep it unchanged. At the end 
one has the number of hard-dimers for a particular configuration 

For fixed N, let Sn be the set of all £jv an d let £n(D) denote the num- 
ber of hard-dimers for £tv- This number is obtained following the steps 1- 
2 above. Then let X = (^n(D))^ nG s n be the sample obtained from this 
numeric counting procedure. The histograms for the standardized sample 
Y = {X — mean (X))/st. deviation (X) are depicted in Figl3]and Fig|4] below. 



ED : = 



Y, N=10 



Y, N=12 
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o 



o 
o 
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Figure 3: Histograms of Y 



Y, N=14 




2 4 6 



Figure 4: Histogram of Y 



Remark 3.1 The red curves in the histograms ofY above represent the stan- 
dard normal distribution. On the other hand, the green curves are obtained by 
smoothing ofY with a Gaussian kernel having standard deviation 0.1. With 
increasing N the green curves become less oscillating and a tendency towards 
a definite distribution becomes visible. Unfortunately, the algorithm becomes 
very time consuming for N larger 14 and it is not clear whether the distribu- 
tions tend to a standard normal distribution. 

In this second part we generate hard-dimers randomly. Roughly speaking, 
this is achieved by first choosing TV" red or blue sites randomly, each with 
probability one half and then for each such configuration dimers are generated. 
Here again randomness enters in that one dimer and the next one are separated 
by some random choice. The algorithm can be implemented in the following 
steps: 

Step 1. Generate a sequence of red and blue sites. The colours appear with 
probability one half. 

Step 2. For the first site, depending on the colour, we choose randomly the 
base elements with attributes equal or right, see FigJSJ For example, if the 
colour is blue, br or be are picked at random with probability one half. 
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These are the building ele- 
ments from which the ran- 
dom hard-dimers are con- 
structed. 

The acronyms mean: re = 
red-equal, rr = red-right, 
rl = red- left, rm = red- 
mixed, and similarly for 
blue. 

Figure 5: Base elements 

Step 3. Suppose we have generated the configuration up to the k th site 
and want to configure the (k + l) th site. Here we have to distinguish the four 
possible combinations of the two colours. Suppose the colours are both blue. 
Then, given the element at the k th site we have to choose the base element of 
the (k + l) th site according to the following table: 



k th site is blue 


(k + l) th site is blue 


br 


bl 


bm 


bm 


be 


random choice € {be,br} 


bl 


random choice € {be,br} 


In the case of a blue and a red site we take: 


k th site is blue 


(k + l) th site is red 


br 


rm 


bm 


rd 


be 


random choice E {re,rr} 


bl 


random choice 6 {re,rr} 



The random choice happens with probability one half. Similarly for the other 
two cases, taking into account the properties of dimers. 

Step 4- For the last site the elements are chosen from the base elements with 
attributes left and equal only. If this is not possible, then a x is being inserted. 
In this way we indicate that the dimer remains unfinished and will therefore 
not be counted in the last step 5, see Figj6j Concretely, we follow the rules: 



— • — rm 


— • — bm 


— rl 


— bl 


•— rr 


• — br 


• re 


• be 
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penultimate site is blue 


last site is blue 


br 


bl 


bm 


X 


be 


be 


bl 


be 


and 


penultimate site is blue 


last site is red 


br 


X 


bm 


rl 


be 


re 


bl 


re 



Similarly for the other two cases. 
Step 5. Count the number of dimers. 



• • • • • • • j^- 

Figure 6: A hard-dimer being omitted 

Now, let Stf denote a random sample of a certain size \S"^\ = m which 
is obtained by repeating the above procedure, steps 1-4, m times. Again, N 
is given by the number of sites. For rj E let rj(d) denote the number of 
dimers counted in r]. Then in Figdthe histograms are shown for two samples 

X = ( V (d)) rieS r N . 
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X, N=15, m=5000 



X, N=30, m=5000 
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Figure 7: Histograms of X 



Remark 3.2 In the histograms for X above, the red curves have been cho- 
sen to have the same mean and standard deviation as the corresponding his- 
tograms. The very good approximation by the histograms makes a central limit 
property visible. 
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